**FULL TITLE** 

ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 



Exoplanet-Induced Chromospheric Activity: Realistic 
Light Curves from Solar-type Magnetic Fields 



o 



Steven R. Cranmer and Steven H. Saar 
- - - Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, 

^ ; Cambridge, MA 02138, USA 

o : 

■ Abstract. There is growing observational evidence for some kind of interac- 

tion between stars and close-in extrasolar giant planets. Shkolnik et al. reported 
^ , variability in the chromospheric Ca H and K lines of HD 179949 and v And that 

I ■ seemed to be phased with the planet's orbital period, instead of the stellar ro- 

tational period. However, the observations also indicate that the chromospheric 
light curves do not repeat exactly, which may be expected for a planet plowing 
through a variable stellar magnetic field. Using the complex solar magnetic field 
I \ (modeled with the Potential Field Source Surface technique) as a guide, we sim- 

^ . ulate the shapes of light curves that would arise from planet-star interactions 

' that are channeled along magnetic field lines. We also study the orbit-to-orbit 

\ variability of these light curves and how they vary from solar minimum (i.e., a 

l/^ , more or less axisymmetric stretched dipole) to solar maximum (a superposition 

^SJ ' of many higher multipole moments) fields. Considering more complex magnetic 

\ fields introduces new difficulties in the interpretation of observations, but it may 

. also lead to valuable new diagnostics of exoplanet magnetospheres. 

O ■ 

2 ■ !• Introduction 

+3 . 

. Many of the ~200 known extrasolar planets are "Hot Jupiters" (i.e., giant plan- 

' ets with star-planet distances less than about 0.1 AU) that have been identified 

•th , from perturbations in their host star's radial velocities. Other ways in which a 

K> \ close-in extrasolar giant planet (CEGP) can influence its star include tidal dis- 

. tortious and magnetic interactions. These effects can produce enhanced stellar 

■ " " ' flare activity (Rubenstein & Schaefer 2000), chromospheric and coronal emis- 

sion (Cuntz et al. 2000), and magnetospheric radio emission (Zarka et al. 2001). 
Tidal and magnetic interactions have been proposed to explain similar phenom- 
ena observed in RS CVn close binary systems (e.g., Simon et al. 1980; Ferreira &: 
Mendoza-Briceno 2005), stars with possibly unseen companions exhibiting "su- 
perflares" (Schaefer et al. 2000), and even in the aurora of Jupiter where there 
are magnetic connections to the inner Galilean moons (Clarke et al. 2002). We 
refer to Saar et al. (2004) for a recent review of observations (or the lack thereof) 
of various kinds of CEGP-star interactions and possible theoretical explanations 
for these effects. 

In this presentation we focus on the chromospheric (Ca H & K line) en- 
hancements on a star due to magnetic interactions between the stellar field and 
the CEGP's magnetosphere. We used the empirically derived solar magnetic 
field, as a function of the 11-year solar cycle, as a proxy for the field geometry 
and field strength of CEGP host stars. We simulated light curves that represent 
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the sum of the rotational modulation of the star's own H&K line emission and 
the added emission from the small set of flux tubes connected to the planet's 
magnetosphere. 

2. Observed Ca H & K Enhancements 

Several early attempts to identify planet-induced chromospheric emission (re- 
viewed by Saar et al. 2004) were unsuccessful mainly because the observations 
were not designed specifically for this purpose. Shkolnik et al. (2003, 2005) used 
high-resolution (A/AA 10^) Ca II H &: K spectroscopy to identify two cases of 
chromospheric enhancement that are phased with the planet's orbit. HD 179949 
(F8 V, Prot ~ 9 days) has a CEGP with Porb = 3.092 days, M s'mi = 0.98 Mjup, 
and semimajor axis d = 0.045 AU. The Ca K enhancement for HD 179949 oc- 
curs at a phase shift A(f) of almost 0.2 ahead of the sub-planet point reaching 
stellar disk-center. This raised suspicions that the star's magnetic field may be 
swept back into a "Parker spiral" at the orbit of the CEGP, but for a solar-type 
magnetic field and wind this is unlikely. On the other hand, HD 179949 is more 
active than the Sun (e.g., Fx/Fq ^ 20) and of earlier spectral type (Tcfj ~ 6170 
K; Valenti &: Fischer 2005) suggesting both the stellar magnetic structure and 
wind could be quite different from the solar case. Shkolnik et al.'s data from 
2003 do not show the strong orbital modulation of the 2001-2002 data, but there 
is still a weak maximum at ^ ~ 0.8 in agreement with the earlier epochs. Similar 
data for v And (F7 V, Prot ~ 14 days, Porb = 4.6 days) reveals a slightly weaker 
H & K amplitude, but a maximum with phase shift ~ 0.5, i.e., when the 
planet is behind the star (see also Harrington et al. 2006). 

3. Theoretical Explanations? 

The periodicity and strength of the HD 179949 and v And variations indicate 
that magnetic interaction is the most likely mechanism. Saar et al. (2004) found 
that other mechanisms such as tidal distortions or lo-likc inductance would 
have produced even stronger modulations in other star-planet systems for which 
no planet-phased variations have been found (e.g., r Boo). In addition, tidal 
distortions would have caused two peaks in the light curve per orbit, rather 
than just one as is observed. The energy released in the Ca H&K enhancement 
(--^10^^ erg/s) is similar to that released in typical solar flares. Ip et al. (2004) 
modeled CEGP magnetospheres and found that this level of power generation 
can be expected from magnetic reconnection in the most likely geometries. 

Mclvor et al. (2006) predicted the shapes of Ca H&K light curves consistent 
with magnetic field lines connecting a planet and an inclined, stretched dipolar 
stellar field. Phase shifts similar to that of HD 179949 can be reproduced by 
misalignment between magnetic and rotation axes. Also, Preusse et al. (2006) 
were able to model similar phase shifts by assuming the CEGP to be inside 
the star's Alfven radius (i.e., where the Alfven speed exceeds the outflow speed; 
r ^ 10 Rq for the solar wind) and that the perturbations travel from the planet 
to the star along the longitudinally varying inward Alfven characteristics. 

Magnetic interaction models depend on the product of the relative velocity 
between the two fields (known from the orbital dynamics), the stellar magnetic 
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field strength (which can be deduced from Zeeman sphtting or estimated from 
known activity scahngs), and the planetary magnetic field strength, which is not 
yet observed, though modeling progresses (e.g., Sanchez-Lavega 2004; Stevens 
2005). 

4. Exploratory Models 

Our goal has been to produce more realistic simulations of Ca K light curves 
for a rotating solar-type star undergoing a CEGP magnetic interaction. We use 
modeled three-dimensional solar magnetic field configurations (as a function of 
the 11-year solar cycle) to simulate both the stellar and planetary light curve 
components. Observed maps of the photospheric magnetic flux were extrapo- 
lated into the extended corona using the Potential Field Source Surface (PFSS) 
method. The corona is assumed to remain current-free out to a spherical "source 
surface" at Rss = 2.5 Rq where the nonradial field components are set to zero, 
simulating the magnetohydrodynamic (MHD) expansion of the solar wind (e.g., 
Schatten et al. 1969; Altschuler & Newkirk 1969; Hoeksema & Scherrer 1986; 
Wang & Sheeley 1990). 

The highest multipole components fall off the most rapidly with increasing 
height, resulting in the photospheric field being much more complex than the 
field at r = Rss- Data and reconstruction coefficients were obtained for 11 
solar rotations (one per year over a solar cycle) from Wilcox Solar Observatory^ 
Somewhat arbitrarily, we used the solar magnetic field from the month of August 
in each of the years of solar cycle 22 (1986-1996); i.e., the specific Carrington 
Rotations 1778, 1792, 1805, 1819, 1832, 1845, 1859, 1872, 1886, 1899, and 1913. 

The following subsections describe how the stellar (§ 4.1) and planetary 
(§ 4.2) light curve components were computed. 

4.1. Stellar Light Curve 

Figure 1 shows the variation of sunspot numbei0 and Ca K-line emission over 
solar cycle 22. The larger number of sunspots at solar maximum corresponds 
directly to the larger number of chromospheric plage and active regions which 
enhance the disk-integrated K-line emission. The Ca K-line data came from 
daily observations at the Sacramento Peak Observatory of the U.S. Air Force 
Phillips Laboratorylfl The Sac Peak K3 index was divided by a fiducial wing-to- 
continuum ratio of 0.40 in order to obtain the same kind of core-to-wing intensity 
ratio used by Shkolnik et al. (2003, 2005). These values were then calibrated 
against the hemisphere-averaged magnetic flux densities (B) (measured at the 
solar surface) over the modeled epochs. We performed a similar fit between 
these two quantities as did Schrijver et al. (1989), and found 

(/corc/Zwing) = 0.13 + 0.024(S)0-63 (1) 



^ http : //sun . Stanford .edu/~wso/ 
' http://sidc.oma.be/sunspot-data7] 

' http://www.ngdc.noaa.gov/stp/SOLAR/ftpcalcium.htnil 
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Figure 1. Top: monthly raw sunspot number over cycle 22. Bottom: Ca II 
K-line ratio of core to wing intensity (see text). Red points in both panels 
define the 11 specific solar rotations that we chose to model in detail. 



where {B), measured in Gauss, was computed from the PFSS models over the 
11 modeled rotations at times corresponding to the K-line observations. This 
relation was then used to compute the theoretical stellar-rotation light curves 
in Figures 3 and 4 from the "exactly known" photospheric magnetic fields. 

4.2. Planetary Light Curve 

To model the enhancement due to planet-star coupling, we assumed the power 
released in magnetic interaction scales as 

^ ~ ^(^^mag)^rel (2) 

where the stellar magnetic field (i?*) and planetary field (-Bp) are measured at 
the planet's magnetosphere, which has a projected area vrr^^g intercepted by 
the stellar field (e.g., Saar et al. 2004). Later, we use rmag ~ 5i?jup, which was 
also assumed by Ip et al. (2004). The two magnetic fields move through one 
another at relative velocity Vrei (~100 km/s for HD 179949). We assumed that 
-Bp, Viaag, and T^ei remain fixed over multiple rotations and the activity cycle, 
and thus that time variations in P are sensitive only to 5*. 

For each of the 11 epochs, we created a fine grid of planetary longitudes 
(300 points around the circular orbit) defined relative to the fixed Carrington 
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Figure 2. Cartoon illustrating the coarse array of magnetic flux tubes traced 
down from the CEGP magnetosphere to the stellar surface. Not to scale. 

longitude coordinates of the modeled magnetic field. For each planetary position, 
we traced an envelope of 5 flux tubes from the planet's orbit (d > Rgs) down 
to the stellar surface. Figure 2 shows the cross-like configuration of flux tubes 
that we used to resolve the magnetosphere; note that the pattern of flux tube 
"footpoints" on the stellar surface is likely to have a more complicated shape. 
We used more than one flux tube so that "bifurcation" events (when the finite- 
sized magnetosphere connects down to two widely separated areas on the stellar 
surface) could be at least approximately resolved. Future work should of course 
use a more finely sampled envelope of flux tubes so that the detailed shapes of 
the resulting star spots can be computed. 

The mapped flux tubes were then input into a time-resolved model of both 
the planet's orbit and the star's rotation. The instantaneous angular positions of 
the star, planet, and observer were tracked, and the angle 9 (measured between 
the stellar surface normal at the flux-tube footpoint and the observer's line of 
sight) was computed as a function of time. The relative contribution from each 
flux tube was thus assumed to be proportional to the stellar fleld strength at 
the magnetosphere (S*) modulated by the spot's projection factor; i.e., 

(^corcZ-^wing) Of ^ (^* COS ^) (summed over 5 footpoints) (3) 

The above procedure contains the assumption that the magnetic "communica- 
tion" between star and planet is essentially instantaneous. This would be a 
valid assumption if the energy were carried by, e.g., suprathermal particles as in 
strong solar flares, but not if the energy had to propagate along more slowly as 
an MHD wave (see Preusse et al. 2006). 

We adopted dynamical parameters roughly analogous to HD 179949 (e.g., 
Pfot = 2.9Porb) not an exact multiple), and we set the overall normalization 
for /coroZ-^wing to resemble the observed amplitude (Shkolnik et al. 2005). This 
arbitrary normalization is of course a weak point of the present models, but 
without a specific predictive model of the magnetic interaction not much more 
was possible. 
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5. Results 

Figures 3 and 4 show the magnetic field configurations and fight curves for the 
11 modeled epochs. The overall polarity "flip" is evident over the 11-year solar 
cycle. Both the complete light curve (star plus planet) and just the stellar rota- 
tional component are shown in order to clearly see the planetary enhancement. 
The light curve is shown in detail over five planetary orbits. The planetary en- 
hancements generally occur when the planet is in front of the star (0 ~ 0) but 
occasional planetary enhancements occur at other phases, or are absent at phase 
zero. Rapid changes sometimes occur over less than 0.05 in phase (i.e., over just 
several hours in time). We also show a much coarser kind of "light curve" that 
was randomly sampled with 200 points spread over a time-span of 200 planetary 
orbits. In these plots, the stellar rotation light curve is not resolved except as a 
fuzzy "band" which the planetary enhancement tends to exceed around phase 
zero. 

Because the amplitudes of both the stellar-rotation light curve and the plan- 
etary enhancements are larger at solar maximum than at solar minimum, there 
is no systematic trend for the planetary enhancement to "poke its head above 
the fuzz" more at either maximum or minimum. Thus, there does not seem to 
be any preferred cycle-phase for the planetary enhancement to be more or less 
distinguishable from the stellar rotational modulation. 

Figure 5 shows some statistical properties of the light curves, with one point 
in each panel characterizing the variations over each of the 11 epochs. The min- 
imum, maximum, and average intensity ratios (/core/-^wing) were determined for 
the separate stellar and planetary light curves and correlated with one another. 
The tendency for the amplitudes of the stellar and planetary light curves to rise 
and fall together was discussed above and is evident in the direct correlations. 
Note specifically the nearly linear relationship between the maximum planet 
enhancement and the total (max — min) amplitude of the stellar light curve. 

6. Conclusions 

Some of the more interesting features of the modeled light curves are summarized 
here: 

1. Because of the complex nature of the multipole fields, the modeled light 
curves do not repeat exactly from orbit to orbit, and sometimes the plan- 
etary enhancement seems to disappear altogether. This may be a possible 
explanation for the 2003 disappearance of the strong orbital modulation 
seen in 2001-2002 for HD 179949 (Shkolnik et al. 2005). 

2. The planetary enhancements arc often non-monotonic in shape (i.e., mul- 
tiply peaked). This occurs because the "spots" on the stellar surface that 
connect to the planet's magnetosphere may be bifurcated or substantially 
swept around in longitude and/or latitude. 

3. For sparsely sampled data, the inferred phase shift between light-curve 
maximum and planetary meridian-passage may range between —0.2 and 
+0.2. 
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Figure 3. For years 1986-1991, the top panel shows the star's magnetic field 
out to i?ss, with positive [negative] polarity in blue [red]. The middle panel 
shows a representative light curve over 5 planetary orbits, and the bottom 
panel shows a coarser random sampling over 200 orbits all phased together. 
In both, the full star+planet light curve is in black (solid curves, filled circles), 
and the stellar rotational component only is in green (dotted curves, X's). 
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Figure 4. For years 1992-1996, the top panel shows the star's magnetic field 
out to i?ss, with positive [negative] polarity in blue [red]. The middle panel 
shows a representative light curve over 5 planetary orbits, and the bottom 
panel shows a coarser random sampling over 200 orbits all phased together. 
In both, the full star+planet light curve is in black (solid curves, filled circles), 
and the stellar rotational component only is in green (dotted curves, X's). 
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Figure 5. Correlations between statistical properties of the stellar (abscissa) 
and planetary (ordinate) light curves over the 11 modeled epochs. 



Better observational phase coverage (even using multiple observatories to obtain 
round-the-clock time series) over consecutive orbits and stellar rotations would 
yield much better tests of the magnetic interaction paradigm. More work needs 
to be done, of course, to model how the stellar magnetic fields differ from the 
fiducial solar case modeled here. Additional kinds of observations, such as X- 
rays that may be phased with the planet's orbit (Saar et al., these proceedings), 
will help constrain the dynamics and energetics of the interaction. 
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